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Abstract 



We perform a detailed study of the decays of the heavier top and bottom squarks 
(t2 and 62) hi the Minimal Supersymmetric Standard Model (MSSM). We show 
that the decays into Higgs or gauge bosons, i.e. ii — > 1 1 + H , A or 
*2 — ► &i,2 + or W + ), and the analogous bi decays, can be dominant in a wide 
range of the model parameters due to the large Yukawa couplings and mixings of t 
and b. Compared to the conventional decays into fermions, such as ti — > t+(xi or <?) 
and t2 —>■ b + Xj~ , these bosonic decay modes can have significantly different decay 
structures and distributions. This could have an important impact on the search 
for ti and 62 and the determination of the MSSM parameters at future colliders. 



The search for supersymmetric (SUSY) particles is one of the most important subjects 
at present and future collider experiments. Future colliders, such as the Large Hadron 
Collider (LHC), the upgraded Tevatron, e + e~ linear colliders, and fi + fi~ colliders will 
extend the discovery potential for SUSY particles to the TeV mass range and allow for a 
precise determination of the SUSY parameters. 

Many phenomenological and experimental studies have been performed for squark (q) 
search In most studies, production and decay of the squarks are studied assuming 
that they decay into fermions, i. e. a quark plus a neutralino chargino (xf), or 

gluino (<?): 

qi^q {,) + (xlxt,ov~g), (1) 

with i,j = 1,2 and k = 1,...,4. However, in the Minimal Supersymmetric Standard 
Model (MSSM) |2| the heavier squarks of the 3rd generation (i.e. stops (t 2 ) and sbottoms 
(b 2 )) can also decay into bosons, i. e. a lighter squark plus a gauge or Higgs boson || [|: 



t 2 -> it Z°, b 2 ^ 6i z° 

t 2 -> &< w+ , b 2 ^ u w- 



and 



(2) 

i 2 - h (h°, H°, A ) , h-^ h (h°, H°, A ) , 



t 2 ^hH + } b 2 -^ U H~ 



(3) 



Here q x (q 2 ) is the lighter (heavier) squark mass eigenstate. The t 2 decays into t\ plus a 
neutral boson in Eqs. (0) and (|3|) are possible in case the difference between m 4 ~ L and mi R 
and/or the t^—tn mixing are large enough to make the necessary mass splitting between 
t\ and t 2 . The same holds for the decays b 2 ^ bi + (Z°, h°, H°, A ). 

In the present article we make a more general analysis than 0, M. We point out that 
the q 2 decays into gauge or Higgs bosons of Eqs. (||) and @ can be dominant in a large 
region of the MSSM parameter space due to large top and bottom Yukawa couplings and 
large t and b mixing parameters. This dominance of the Higgs/gauge boson modes over 
the conventional fermionic modes of Eq. ([]]) could have a crucial impact on searches for 
t 2 and b 2 at future colliders. 

First we summarize the MSSM parameters in our analysis. In the MSSM the squark 
sector is specified by the squark mass matrix in the basis (q L , q R ) with q = t or b |j| § 

Ml = ( < ^ ) (4) 



a q m q m~ R 



with 



m\ L = Ml + m 2 z cos2(3(II L -e q sm 2 e w ) + m 2 q , (5) 

m 2 g R = M 2 {i] p } + m% cos 2(3 e q sin 2 9 W + m 2 q , (6) 

j (A t - fx cot p) m t (q = t) 

aqmq " 1 (A-/itan/3)m b (q = b). UJ 
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Here If is the third component of the weak isospin and e q the electric charge of the 
quark q. Mq^jj and A tjb are soft SUSY-breaking parameters, /i is the higgsino mass 
parameter, and tan/3 = v 2 /vi with v\ [v 2 ] being the vacuum expectation value of the 
Higgs H\ {H®). From renormalization group equations we expect that Mq, Mq, and 
of the third generation are different from those of the first two generations and from each 
other. Diagonalizing the matrix (||) one gets the mass eigenstates q x = q L cos9q + q R sin9q, 
q 2 = — ?i,sin0q + QR cos @q with the masses m^, rriq 2 (m^ < m ? - 2 ) and the mixing angle 
9g. As can be seen, sizable mixing effects can be expected in the stop sector due to the 
large top quark mass. Likewise, b^-b R mixing may be important for large tan f3. 
The properties of the charginos xf (i = 1,2; m~± < Tn~±) and neutralinos x°k {k = 1, ...,4; 
m^o < ... < rn^o) are determined by the parameters M, M' , fi and tan/5, where M and M' 
are the SU(2) and U(l) gaugino masses, respectively. Assuming gaugino mass unification 
we take M' = (5/3) tan 2 6 W M and rrig = (a s (m g ) / 'a 2 )M with rrig being the gluino mass. 
The masses and couplings of the Higgs bosons h°, H°, A and H^, including leading 
Yukawa corrections, are fixed by rriA, tan/3, fi, m t , m b , Mq, Mjy, Mg, A t , and A b . 
H° (h°) and A are the heavier (lighter) CP-even and CP-odd neutral Higgs bosons, 
respectively. For the Yukawa corrections to the h° and H° masses and their mixing angle 
a we use the formulae of Ref. J7J. For H ± we take m 2 H± = m\ + m^. 

The widths of the squark decays into Higgs and gauge bosons are given by = 1,2; 
k = l..A) i: 

Hi v Hi 

Here H k = {h°, H°, A , H^} and V = {Z°,l¥ ± }. The G ijk denote the squark couplings 
to Higgs bosons and c^y those to gauge bosons. k^ x = «(m| , m^, nix) is the usual 

kinematic factor, n(x,y, z) = (x 2 + y 2 + z 2 — 2xy — 2xz — 2yz) 1 / 2 . (Notice that r(^ — * 
qf 1 H k ) is proportional to k whereas r(g j — >■ V) is proportional to k 3 .) The complete 
expressions for and well as the widths of the fermionic modes, are given in 

0, |H . The leading terms of the squark couplings to Higgs and vector bosons are given in 
Table 1. The Yukawa couplings ht,b are given as 

h t = g m t /(V2m w sin /3), h b = g m b /(y/2m w cos (3) . (9) 

As can be seen, the q x q 2 Z® couplings take their maximum for full q^-q^ mixing [9q — > 7r/4 
or 37r/4) and vanish in case of no mixing. The reason is that the Z° couples only to q^q^ 

and q R q R - On the other hand, the W ± couples only to the left components of the squarks. 

(') 

In contrast to that, Higgs bosons couple mainly to q L -q R combinations. These couplings 
are proportional to the Yukawa couplings h t)b and the parameters A tyb and /i, as can 
be seen in Table 1. Notice here that the qiq 2 h° and q^H couplings have a factor 
cos 29g (which decreases with increase of the g-mixing) while the qiq 2 A° couplings do not 
depend explicitly on the squark mixing angles. Hence the ^1^2^.° coupling (and the bib 2 A° 
coupling for large tan/3) can be especially strong in case A t (A b ) and fi are large. Notice 
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t x t 2 Z Q ~ g sin26»i 

tit 2 h° ~ ht (fM sin a + At cos a) cos 2% 

tit 2 H° ~ (/i cos a — v4 4 sin a) cos 26> 4 - 

tit 2 A° ~ (/i sin /? + A cos /?) 



6i6 2 Z° ~ 5 sin 2^ 

b x b 2 h° ~ /ib (/i cos a + Afe sin a) cos 2^ 
bib 2 H° ~ /i b (/isina — v4 b cosa) cos 2^ 



&i& 2 A° ~ /i b (/i cos /3 + A b sin /?) 



t 2 b 2 H 



-h t (// sin /3 + 

3 y \ -sm^-cos^ 



h t (/isin/3 + A t cos/5) sin 6*^ cos ^ + ^ (/xcos/? + A^sm/?) cos^sin^ 
t\b 2 H ± ~ — /i t (/isin/? + A t cos/3) sin^ sin^ + /if, (/icos/3 + Afosin/?) cos^cos^ 
/i t (/i sin j3 + A t cos /3) cos 6^ cos 6^ — /i;, (yU cos (3 + A b sin f3) sin 6*^ sin 0^ 
> — /^t (/i sin (3 + A t cos /3) cos % sin ^ — h b (fi cos /3 + A b sin /?) sin 6$ cos 



— cos 6i sin ^ 

shift- sin #r / 

1 b ' y 



Table 1: Squark couplings to Higgs and vector bosons (leading terms). 



also that the squark mixing angles themselves depend on A tt b, At, and tan/3. Moreover, 
Ql~Qr mixing enhances the splitting of the q mass eigenvalues, which in turn can have an 
important influence on the phase space of the q 2 decays. 

We now turn to the numerical analysis of the t 2 and b 2 decay branching ratios. For this, 
we calculate the widths of all possibly important 2-body decay modes (3-body decays || 
are negligible in this study): 



*2 -> tg, txt bxj, t x Z\ hW + , t x h\ t!# u , t x A u , hH 



b2 - b~g, bxl txj, b i Z °> tiW~, hh°, hH , b x A°, UH- 



We take m t = 175 GeV, m b = 5 GeV, m z = 91.2 GeV, sin 2 6 W = 0.23, m w = m z cos 9 W , 
a{m z ) = 1/129, and a s (m z ) = 0.12 [with a s (Q) = 127r/((33-2n/) ln(Q 2 /A^ )), n f being 
the number of quark flavors]. In order not to vary many parameters we choose Mq = 
|Mj) = j^Mp and A t = A b = A for simplicity. Moreover, we fix M = 300 GeV (i.e. 
rrig = 820 GeV) and = 150 GeV. Thus we have as free parameters Ma, /i, tan/3, and 
A. In the plots we impose the following conditions: 

(i) m-± > 100 GeV, 

(ii) m^o > 70 GeV, 

(iii) m hh > m x?' 

(iv) m h o > 80 GeV, 

(v) Ap (t — b) < 0.0016 || using the formula of flQ|| , and 
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(vi) A\ < 3 (M~ + M? + m 2 ^) and A 2 < 3 (M| + M? + m^) with 

m^ 2 = (m\ + m|) cos 2 /3 — | m 2 z and = (m 2 ^ + m|) sin 2 /3 — ~ 
(approximately necessary conditions to avoid colour and electric charge breaking 
global minimum ]TT||). 

Conditions (i)-(iv), along with = 820 GeV, satisfy the experimental bounds on xt, xl, 
t, b, and h° from LEP2 [|i~2| , [T3| , |14"|| and Tevatron jnj. Conditions (v) and (vi) constrain 
the t and b mixings significantly. 

In Fig. |l] we plot the contour lines for the branching ratios of the Higgs boson modes and 
the gauge boson modes, i. e. BR(g 2 q^H) = £ BR[ q 2 q x + (h°, H°, A ), cf 12 + ] 
and BR(g 2 -> q^V) = £BR[ q 2 -* q x + Z°, q' l2 + W ± ] with g = t or 6, in the fi-A 
plane for Mg = 500 GeV and tan/3 = 3. We see that the i 2 and 62 decays into bosons 
are dominant in a large region of the MSSM parameter space, in particular for large |/i| 
and/or \A\. Note here the dependence on the signs of A and /1. We have obtained a 
similar result for large tan f3. 

In Fig. H we show the individual branching ratios of the t 2 and 62 decays into bosons as 
a function of // for Mg = 500 GeV, A = 600 GeV, (a) tan/3 = 3 and (b,c) tan/3 = 30. 
(We plot only branching ratios larger than 1%.) In case of tan/3 = 3 (Fig. 2a) the 
bosonic decays of t 2 are dominant (BR > 50%) for fi < — 400 GeV because (i) (At — 
/icot/3)m t , the off-diagonal element of the stop mass matrix, is large enough to induce 
the necessary mixing and mass splitting for the stops, and (ii) for relatively large the 
decays into higgsino-like neutralinos (x° 4 ) and chargino (xf) are kinematically suppressed 
or forbidden. The branching ratio of the tiZ° mode has its maximum at fi = —600 GeV 
where = 332 GeV, m 4 ~ 2 = 628 GeV, and 61 ~ 130°. For further decreasing /i the decay 
into t\A° quickly gains importance as the tit 2 A° coupling is ~ h t (0.3A t + 0.95 /i). The 
b\W + mode has a branching ratio of 10% to 16% for fi < — 400 GeV because 61 ~ bi 
and m~ b ~ 500 GeV. For positive /1 the bosonic t 2 decays are kinematically suppressed 
or even forbidden due to insufficient stop mass splitting. Similar arguments hold for the 
b 2 decays. Here only the decays into t\H~ and t\W~ are kinematically accessible, which 
can be seen in Figs, lc and Id. 

For tan (3 = 30 (Figs. 2b and 2c) there is a large mixing also in the sbottom sector (which is 
proportional to /ztan/3 while % ~ 128°, see Eq. (|7|)). Hence the t 2 and b 2 decay branching 
ratios show a much weaker dependence on the sign of /1. For > 400 GeV the decays 
into bosons clearly dominate. Notice the importance of the decays t 2 — * bi + (H + , W + ) 
in Fig. 2b and b 2 — ► b\ + (A , Z°) in Fig. 2c. This is mainly due to the large mixing (and 
hence large mass splitting) in the sbottom sector. Notice also that, in general, when the 
decays into Z° and A are kinematically allowed, those into h° and H° are also possible. 
The latter decays are, however, practically negligible in this example because they are 
suppressed by a factor cos 2 26 q. (Here note that BR(t 2 t\h°) and BR(t 2 — * ii-ff ) can 
be about ~ 10% for other values of the MSSM parameters; see e.g. Fig. 4.) 
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In Fig. |3] we show the tan f3 dependence of the bosonic £2 and 62 decay branching ratios 
for Mq = 500 GeV, A = 600 GeV, and \i = —700 GeV. Here it can be seen explicitly how 

the q 2 decays into the lighter sbottom b\ plus a gauge or Higgs boson become important 
with increasing tan/?. 

In Fig. H we show the Mq dependence of the branching ratios of the bosonic t% and 
b 2 decays for A = 400 GeV, /i = —1000 GeV, and tan/5 = 3. In this case we have 
(m^,m-+,mj) = (150, 302, 820) GeV. We see that the bosonic modes dominate the t 2 

and b 2 decays in a wide range of Mq. (Notice that the decay into a gluino is dominant 
above its threshold.) 

For large rriA the decays into H°, A , and H are phase-space suppressed. However, 
the remaining gauge boson modes can still be dominant. For Mq = 500 GeV, fi = 1000 
GeV, A = 1000 GeV, and tan/3 = 3 we have, for instance, (m^, m 4 ~ 2 , m~ bi , m~ b2 , m~+, 

m -o) = (365, 610, 501, 558, 294, 148) GeV with BR(£ 2 -> qH) = (68, 58, 32, 1)% 
and BR(t 2 -> qV) = (18, 24, 39, 57)% for m A = (110, 200, 240, >250) GeV. As for the 
dependence on the parameter M, our results do not change significantly for smaller values 
of M, as long as decays into gluino are kinematically forbidden. 

Let us now discuss the signatures of the t 2 and b 2 decays. Typical signals of the decays 
into bosons (Eqs. (|2|) and (|3|)) and of those into fermions (Eq. ([!])) are shown in Table 1 
of Ref. |16[]. In principle, the final states of both types of decays can be identical. For 
example, the final state of the decay chain (a) t 2 — > t\ + (h°, H°, A or Z°) — > (txi) + 
(bb) — > ibqq'xi) + (bb) has the same event topology as that of (b) i 2 —* t + x 2 3 4 — > t 
+ ((h°,H°,A° or Z°) + xi) -> t + (bbxl) -> ibqq') + (66x2)- Likewise (c) t 2 ^ ~b 1>2 
+ (H + or W + ) — * (bxl) + (qq') has the same event topology as (d) t 2 — > b + xl,2 ~^ 
b + ((H + or H^ + ) + Xi) — > b + (qq'xi)- However, the decay structures and kinematics 
of the two modes (a) and (b) ((c) and (d)) are quite different from each other, since 
the Xi (supposed to be the lightest supersymmetric particle (LSP) and hence a missing 
particle in case of R-parity conservation) is emitted from t\ and X234 (^1,2 an d Xi^)^ 
respectively. This could result in significantly different event distributions (e.g. missing 
energy-momentum distribution) of the q 2 decays into gauge or Higgs bosons compared to 
the decays into fermions. Hence the possible dominance of the former decay modes could 
have an important impact on the search for t 2 and 62, and on the measurement of the 
MSSM parameters. Therefore, the effects of the bosonic decays should be included in the 
Monte Carlo studies of i 2 and b 2 decays. 

In conclusion, we have shown that the t 2 and b 2 decays into Higgs or gauge bosons 
(such as t 2 — >• t\ + (h°, H°, A or Z°) and t 2 — > b± .2 + {H + or W + )) can be dominant in a 
wide range of the MSSM parameter space due to large Yukawa couplings and mixings of 
t and b. Compared to the conventional fermionic modes (such as t 2 —* t + (x® or g) and 
t 2 — >• b + Xj)i these bosonic decay modes can have significantly different decay structures 
and distributions. This could have an important impact on the searches for £2 and &2 and 
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on the determination of the MSSM parameters at future colliders. 
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Figure 1: Branching ratios (in %) of t 2 and b 2 decays in the \i-A plane for A t = Ab = A, 
M Q = 500 GeV, Mjj = AAA GeV, M b = 556 GeV, M = 300 GeV, m A = 150 GeV, and 
tan/3 = 3. (a)EBR[f 2 ^f 1 +(/i° L //°,4 ), &i, 2 +if+ ], (b) E BR[t 2 -> ti + Z°, b lt2 +W + ], 
(c) BR[ 6 2 — ► t\ + -ff - ], (d) BR[ b 2 — > £1 + ]. In the dark grey areas the decays are 
kinematically not allowed; the light grey areas are excluded by the conditions (i) to (vi) 
given in the text. 
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Figure 2: fi dependence of ti and 62 decay branching ratios for Mq = 500 GeV, My = 
AAA GeV, M b = 556 GeV, A t = A b = 600 GeV, M = 300 GeV, m A = 150 GeV, (a) 
tan/5 = 3, (b) and (c) tan/3 = 30. The grey areas are excluded by the conditions (i) to 
(vi) given in the text. 
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Figure 3: tan/5 dependence of t 2 and 6 2 decay branching ratios for Mq = 500 GeV, 
My = 444 GeV, M b = 556 GeV, A t = A b = 600 GeV, /i = -700 GeV, M = 300 GeV, 
and rriA = 150 GeV. 



11 




Figure 4: Mq dependence of t 2 and b 2 decay branching ratios for Mjy = | Mq, = 
fMg, A t = A b = 400 GeV, // = -1000 GeV, M = 300 GeV, m A = 150 GeV, and 
tan p = 3. 
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